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The physical consequences of the spin-triplet, even-parity pairing that has been predicted to exist in 
disordered two-dimensional electron systems are considered in detail. We show that the presence of 
an attractive interaction in the particle-particle spin-triplet channel leads to an instability of the nor- 
mal metal that competes with the localizing effects of the disorder. The instability is characterized 
by a diverging length scale, and has all of the characteristics of a continuous phase transition. The 
transition and the properties of the ordered phase are studied in mean-field theory, and by taking 
into account Gaussian fluctuations. We find that the ordered phase is indeed a superconductor with 
an ordinary Meissner effect and a free energy that is lower than that of the normal metal. Various 
technical points that have given rise to confusion in connection with this and other manifestations 
of odd-gap superconductivity are also discussed. 



I. INTRODUCTION 

The possibility of superconductivity with a gap func- 
tion that is an odd function of time or frequency has been 
the subject of some attention lately. The concept was in- 
troduced by Berezinskiity in the context of 3 He, and it was 
more recently revived in connection with-two-dimensional 
(2-D) electron nuidsJn semiconductorsja and with high- 
T c superconductors.0 In Berezinskii's original work, as 
well as in Refs. 0,f|, the case of spin-triplet, even-parity 
pairing was considered, while Refs. || also discussed the 
case of spin-singlet, odd-parity pairing. In either case, 
the gap function being an odd function of the frequency 
ensures that the Pauli principle is obeyed. 

In these references, some properties of odd-gap super- 
conductors have been explored, but no complete analysis 
of the phase transition, or of the superconducting phase, 
has been given. Furthermore, several fundamental ques- 
tions concerning the stability of such a superconducting 
phase, and whether an odd-gap superconductor can in- 
deed be superconducting, have led to considerable con- 
fusion. 

In the present paper we discuss these issues. For defi- 
niteness, we will analyze the mechanism for spin-triplet, 
even-parity superconductivity discovered in Refs. |], but 
many of our conclusions apply to other odd-gap super- 
conductors as well. The paper is organized as follows. 
In Sec. |l| we give a field-theoretic formulation of the 
problem and discuss the stability properties of the field 
theory. In Sec. Ill we show that in the presence of an 



attractive interaction in the particle-particle spin-triplet 
channel, there is an instability of the normal conduct- 
ing phase. This instability has all of the properties of 



a second order phase transition, with a diverging length 
scale, a diverging order parameter susceptibility, etc. In 
Sec. |[y] we develop and analyze a mean- field theory of 
the transition. In particular, we determine the critical 
behavior on either side of the transition. We then deter- 
mine the properties of the ordered phase in a Gaussian 
approximation. We show that the ordered phase is a su- 
perconductor, with a Meissner effect, and a real part of 
the conductivity that has a delta-function contribution, 
just as in BCS superconductors. In Sec. [v| we conclude 
with a summary and a general discussion of our results. 



II. FIELD-THEORETIC FORMULATION OF THE 
PROBLEM 

A. Q-matrix theory for fermions 

Our starting point is a general field theory for elec- 
trons. For anv, fcrmionic system, the partition function 
can be writtcnla 



Z = J D$,1>] exp(ShM) 



(2.1a) 



where S is the action in terms of the fcrmionic (i.e., 
Grassmann valued) fields -0 and i/j. We consider an ac- 
tion that consists of a free-fermion part So, a part S<ji s 
describing the interaction of the electrons with quenched 
disorder, and a part Si n t describing the electron-electron 
interaction, 



S — Sq + S, 



dis 



Sir 



(2.1b) 
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Each field -0 or t/5 carries a Matsubara frequency index 
n, a spin index a =T, I, and, if the quenched disorder is 
dealt with by means of the replica trick, a replica index 
a. For our purposes it is useful to introduce a matrix of 
bilinear products of the fermion fields, 



/ 



B\2 = 



-0) 1T -02T 



-i>w>p2i 



Hi 



V'itV^t J 

(2.2) 



where all fields are understood to be taken at position x, 
and 1 = (ni,ai), etc. The matrix elements of B com- 
mute with one another, and are therefore isomorphic tp, 
classical or number- valued fields that we denote by QE 
This isomorphism maps the adjoint operation on prod- 
ucts of fermion fields, which is denoted above by an over- 
bar, on the complex conjugation of the classical fields. 
We use the isomorphism to constrain B to the classical 
field Q, and exactly rewrite the partition function] 

Z = J e s ^l J D[Q]5[Q-B] 

' D®,1>] e s ^ { D[Q]D[1] e^'IMQ-s)] 



(2.3) 



= / D[Q]D[A] 



A[Q,A] 



A is an auxiliary bosonic matrix field that serves to en- 
force the functional delta-constraint in the first line of 
Eq. (2.3), and the last line defines the action A. The 



matrix elements of both Q and A are spin-quaternions 
(i.e., eleme nts of Q x Q with Q the quaternion field). 
From Eq. (2.2) we see that expectation values of the 
Q matrix elements yield local Green functions, and Q- 
Q correlation functions describe four-fermion correlation 
functions. The physical meaning of A is that its expecta- 
tion value plays the role of a self energy (see Ref. and 
Sec. below). 

It is convenient to expand the 4x4 matrix in Eq. (|2. 
in a spin-quaternion basis, 



Qi2(x) = ^2 ( 7V 



&ia(x) 



(2.4) 



r\i=0,3 



and analogously for A. Here to = So = I2 is the 2x2 unit 
matrix, and tj = — Sj = —i<Jj, {j = 1,2,3), with CTi.2,3 
the Pauli matrices. In this basis, i = and i — 1,2,3 de- 
scribe the spin singlet and the spin triplet, respectively. 
An explicit calculation reveals that r — 0, 3 corresponds 
to the particle- hole channel (i.e., products ipip), while 
r — 1,2 describes the particle-particle channel (i.e., prod- 
ucts ipip or il>ip)- We will be particularly interested in the 
matrix elements of \Q, for which the isomorphism ex- 
pressed in Eq. (O) readaS 



1Q12W - g [V'itMV^tM -^u( x )^2i(x) 

+ u (x)02i(x) - & T (x)& t (x)] . (2.5) 

From the structure of Eq. ([^) one obtains the,-following 
formal symmetry properties of the Q matrices, 



(2.6a) 
(2.6b) 
(2.6c) 
(2.6d) 
(2.6e) 
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rQ\2 
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Here the star in Eq. ( 2.6ej ) denotes complex conjugation. 

By using the delta constraint in Eq. (2.3) to rewrite 
all terms that are quartic in the fermion field in terms of 
Q, we can achieve an integrand that is bilinear in ip and 
tp. The Grassmannian integral can then be performed 
exactly, and we obtain for the action A 



A[Q, A] = „4 dis + A int + I Tr In ( G 



dx tr ( A(x) Q(x) 



iA 



Here 



Go 1 = -d r + d 2 j2m 



(2.7a) 



(2.7b) 



is the inverse free electron Green operator, with d T and <9 X 
derivatives with respect to imaginary time and position, 
respectively, m is the electron mass, and /i is the chemical 
potential. Tr denotes a trace over all degrees of freedom, 
including the continuous position variable, while tr is 
a trace over all those discrete indices that are not ex- 
plicitly shown. The electron-electron interaction .4i nt is 
conveniently decomposed into four pieces that describe 
the interaction in the particle-hole awl particle-particle 
spin-singlet and spin-triplet channels .0 For the purposes 
of the present paper, we need only the particle-particle 
spin-triplet channel interaction explicitly. In Ref. ^ it 
was shown that in any quenched disordered, interacting 
electron system, there is an attractive interaction in the 
particle-particle spin-triplet channel of the form 



A 1 



p-p,t 

int 



x K ni . n2 . n3trl4 ^2 *^2*^2rQ ni n 2 rQn 3 n 4 • ( 2 -8a) 
r— 1,2 i— 1 a 

In D = 2, the effective interaction potential is 
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K ni ,n2;n 3 ,n4 — . (^n 1 ,n 2 ;™3,™4 ,1 " i •'» i : ■ 



Kn 



^ni ,ri2 ;ri4 .713 K n2 ni - ni _ n3 ) , (2.8b) 



whereB 
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K 



ni.n2;i3,i4 



y In 



n 2 - n 3 



(2.8c) 



with a positive coupling constant y > that depends 
both on the disorder strength and on the coupling con- 
stants in the particle-hole interaction channels. Although 
both the disorder and the particle-hole channel interac- 
tions are necessary to produce the particle-particle spin- 
triplet interaction, none of the points to be investigated 
in this paper qualitatively depends on either one of them 
other than through the existence of K. For the Gaus- 
sian theory that we will consider, the only other effect 
of the disorder that is relevant to our dicussion is that it 
replaces some free-electron correlation functions by dif- 
fusive ones, which changes the exponents in certain scal- 
ing relations. For notational simplicity, and in order to 
keep our discussion technically as simple as possible, in 
what follows we therefore neglect all contributions to the 
action that are nonessential for our p urposes. We thus 
work with a system given by Eq. ( 2.7a ) with _4dis = and 
Ant = -4fnt P ' and drop the replica indices on all fields. 
In cases where diffusive correlations make a difference, we 
will mention this explicitly and restore diffusive scaling. 
Our restriction to a Gaussian approximation purposely 
neglects the localizing effects of th e disorder. We will 



come back to this poin t in Sec. |V C 
Note that Eq, 



(2- 



means that in time space, the in- 
teraction between superconducting fluctuations is long- 
ranged. Consequently, the critical behavior, and in par- 
ticular the critical exponents discussed in Sees. [II and 
IV A below, depend explicity on the detailed form of the 
kernel. However, our qualitative results, in particular 
the spontaneous symmetry breaking and the existence of 
an ordinary quantum critical point that marks the on- 
set of superconducting long-range order, we expect to be 
generic. We also note that in Eqs. (2.8) above, as in Ref. 



m, we have neglected any wavenumber dependences of the 
effective interaction potential. Any nontrivial, i.e. non- 
analytic, wavenumber dependence would correspond to 
an interaction that is long-ranged in real space. Such 
long-ranged interactions are known to stabilize mean- 
field critical behavior Jill and will thus increase the range 
of validity of our mean- field theory. We therefore do not 
expect our neglecting the wavenumber dependence of K 
to qualitatively affect our results. 



B. The Fermi-liquid saddle point, and Gaussian 
fluctuations 

In Ref. it was shown that the above Q-A field the- 
ory possesses a saddle-point solution that describes a free 
Fermi gas (or a disordered Fermi liquid if we had not 
dropped the disorder and particle-hole channel interac- 
tion contributions to the action). That is, the saddle- 
point equations 



SA 
SQ 



Qsp.Asp 



SA 
SA 



= 



Q H p,A Bp 



are solved by the ansatz 

r Q 12 (x) = 8 12 SrO SiQ Q ni , 



;a 12 ( x ) 



sp 



8 12 8 r Q 8iQ A ni , 



with 



A„=0 . 



(2.9) 

(2.10a) 
(2.10b) 

(2.11a) 
(2.11b) 



The saddle-point Green function G° is simply the one for 
free electrons, 



G°(p) = (io; n -ep)- 1 , 



(2.12) 



with £ p = p 2 /2m — /it. 

We next consider the Gaussian fluctuations about this 
saddle point. Proceeding as in Ref. 0, we write 



Q = Qs P + SQ 
A = A sp + SA. 

and introduce a new field A by 



Al2 = - ¥?12 A12 



!12 



(2.13a) 
(2.13b) 



(2.14a) 



with 



^nm(k) = G sp(P)W n )G S p(p + k,cJ m ) , (2.14b) 



a generalized Lindhard function. A has been chosen so 
that Q and A decouple in the Gaussian action. Expand- 
ing to second order in the fluctuating fields SQ and SA, 
we obtain the Gaussian action 



Ac 



EEE 



^(k) [;(*A) ia (k) 



v 

k 1,2 r,» \ + / r 

x *(*A)i 2 (-k) - *(<fQ)i 2 (k)*(aQ) ia (-k)] 



where the symbol 



+Af-^[8Q} , (2.15) 



is equal to +1 for r — 0, 3, and 



-1 for r = 1,2. 

Keeping the disorder part of the actio n resu lts in the 
inverse Ligdhard function (^^(k) in Eq. (2.15) being re- 
placed byEI 



©r 2 1 (k) = ¥T 2 1 (k)-l/ 7 riV P T , 



(2.16a) 
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with r the elastic scattering mean-free time. For small where the matrix M has the structure 



wavenumbers and small u> ni 
diffusive, 



X>i 2 (k) 



ttN v 



Dk 2 



with ni7i2 < 0, T> is 



(2.16b) 



with D the Boltzmann value of the diffusion constant. 
Whenever it makes a difference in the results, we will use 
the diffusive propagator T> instead of the free fermion 
propagator ip. 

In order to completely define the field theory, one also 
must specify the integration contours for the functional 
integral, Eq. (2J3). The most obvious choice would be to 
i nteg rate over the space of all matrices that obey Eqs. 
( |2.6|) . However, t his c re ates t he following problem. Ac- 
cording to Eqs. (2.6d), (2.6e), \Q n _ l _ n is real, and so 



are all other particle-particle spin-triplet components of 
Qn-l,-n- However, according to Eq. ( |2.15[ ), this direc- 
tion in the complex Q-plane is unstable, and the direction 
of steepest descent is along the imaginary axis.Eil In or- 
der to do the integral by the saddle-point method, the 
contour for the integration over the 12 Q must therefore 
be deformed so that it passes through the saddle point in 
the direction of the imaginary axis. At least within per- 
turbation theory, this just amounts to formally doing the 
Gaussian integral without worrying about convergence 
problems. A necessary condition for this procedure to be 
valid is that the resulting theory reproduces the pertur- 
bative results obtained within the underlying fermionic 
theory. An easy and convenient check is provided by the 
so-called weak-localization correction to the conductivity 
of non-interacting electrons. For our purposes, a simple 
structural check suffices and is provided in Appendix |A|. 
The check is affirmative, and suggests that we can safely 
ignore convergence questions relating to the Gaussian in- 
tegrals that occur in our field theory. A related point will 
become important in Sec. IV below. 



III. INSTABILITY OF THE NORMAL METAL 

A. The Gaussian eigenvalue problem in the normal 
phase 

We now turn to the full Gaussian action in an ex pan- 
sion about the Fermi liquid saddle-point, Eq. ( 2.15 ). In- 
tegrating out the auxiliary field A just contributes a mul- 
tiplicative constant to the partition function, and we are 
left with a quadratic form in SQ. Let us consider the 
particle-particle spin-triplet part of the action, e.g., the 
channel r — i — 1. Since the action is diagonal in all 
indices except the frequency, we can drop the r and i 
indices as well as the wavenumber dependence, and con- 
sider a quadratic form 



1,2,3,4 



'12 



M 



12.34 



SQ 



31 



(3.1a) 



M 



12,34 



= -a\ 2 #13 #24 - h- 



3,2-3 01+2,3+4 



(3.1b) 



Here b is proportional to the interacti on po tential K, 
Eq. ( 2.8c ), and ayi is either <fii 2 ' Eq. ( 2.14b ), or , 
Eq. (2.16a), depending on whether or not we keep the 
disorder explicitly. It is useful to consider M as a matrix 
with composite indices (1, 2) and (3, 4), and to study the 
eigenvalue problem 



l',2' 



M 



12,1'2' jvn 



A/i 



(3.2) 



with eigenvalues A and eigenfunctions fyi- For reasons 
that will soon become apparent, it is advantageous to 
transform to a different basis of eigenfunctions g defined 
by 



312 



J2 b 

3,4 



1-3,2-3 "1+2,3+4 /34 



(3.3) 



In terms of the g, the eigenvalue problem reads 
912 -- 



Z) „.. ''' x gi',-i'+(i+2) • (3.4a) 



a l',-l' + (l+2) + ^ 



Now we write uj ri2 = —uj ni + Q n , analytically continue to 
real frequencies at T — 0, and put the wavenumber de- 
pendence back in. Then we obtain the eigenvalue equa- 
tion in the form 



gx(uj;p,n) = -y 



dx In 



1 



LO — X 



a(x; p, fi) + A 



d — x + n 

gx(x;p,n) . (3.4b) 



Here y is the coupling constant of Eq. (2.8c), scaled by 
an appropriate factor. In order to check for an insta- 
bility of the metallic phase, we need to look for a zero 
eigenvalue. From the structure of Eq. ( 3.4b| ) it is clear 
that the first zero eigenvalue appears for p = £1 = 0. We 
therefore specialize to zero momentum and external fre- 
quency. We then have a(x; p = 0, il = 0) oc \x\, irrespec- 
tive of whether we use a — ip^ 1 or a = V~ x . Absorbing 
a constant factor into the eigenvalue A and into the cou- 
pling constant y, we then get the eigenvalue equation in 
the form 



dx 
x + X 



In 



X + LU 



X 



LO 



9x(x) 



(3.5) 



with g x (u>) = g\{uj; p = 0, Q = 0). Notice that Eq. fl3.5| ) 
is a generalization of the gap equation in Ref. |4j: The 
critical eigenfunction, 5^=0(^)1 obeys the same integral 
equation as the critical gap function. This is the first in- 
dication that the long-range order implied by a nonzero 
gap function in Ref. H can indeed be understood in terms 
of a conventional continuous phase transition that is trig- 
gered by an instability of the metallic phase. 
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B. Solution of the eigenvalue problem 



The integral equation, Eq. (3.5), is very hard to solve 
analytically. We therefore make the same approximation 
as was done for the gap equation in Ref. [|, namely re- 
placing the logarithmic kernel by a rational one with a 
similar overall behavior: 



In 



X + LU 



X — Ul 



— e{x~to) + — e{Lu~x) . (3.6) 

Ul ' X 



This allows to rewrite the integral equation as an ordi- 
nary differential equation, 



d ( \ 



1 d 
u) du> 



.9a H 



1 



u){uj + A) lu 2 



g\{u) = o. 

(3.7) 



It is useful to rewrite this ODE in the form 



d 2 d 

-j^ w(z) + p(z) — w(z) + q(z) w(z) = , (3.8a) 



with 



p(z) = - 



t-1 1-t 



1 



(3.8b) 



(3.8c) 



where t = 1 — 4y, and the solution w(z) of the ODE 
determines the eigenfunction via 



g x (uj) = \ 2 w(-ui/\) 



(3.8d) 



The ODE, Eq. ( 3.8a ), is Fuchsian with three regular 
singular points at z = 0, z = 1, and z — oo. ItrGan thus 
be transformed into a hypergeometric equation.113 Before 
we look at the general solution, let us con sider the case 
A = 0, where the two solutions of Eq. (3.7) are 



o (1) 



-Vt 



n {2) 

9x=o - u 



(3.9) 



In order to obtain a solution that is well-behaved every- 
where, we therefore must require t(X = 0) = 0. More 
generally, the requirement of a well-behaved eigenfunc- 
tion leads to t being a functio n of A, as we now proceed 
to show. The solution of Eq. ( 3.8a ) that is well-behaved 
for to — > oo can be written 



-Vt 



1 



xF 1 



%/7,2 + Vt;3; 



1 



1 + A/w 



(3.10) 



with F a hypergeometric function. For later reference 
we also list the large-frequency behavior of the original 
eigenfunction /, related to g by Eq. (3.3). By an analysis 
analogous to the one above for g one finds 



f\{ui — > oo) oc (uj/X) 



-i-Vt 



(3.11) 



C. The correlation length, and the exponents v, 77, 
and 7 



Let us consider the eigenfunction, Eq. ( 3.10| ), in the 
limit A — > 0. Anticipating that t — > as well in that 
limit, we expand in t and find asymptotically 



5a^oM = -(w/A)" 



(3.12) 



We know, however, that in this limit the eigenfunction 
must merge with the solution g^l , Eq. (3.9). This yields 
A as a function of t for asymptotically small t, 



A(t -» 0) oc t 1 ' 2 ^ 1 



(3.13) 



Furthermore, we know that at finite momentum p, both 
the eigenvalue A and p 2 appear additively in the eigen- 
value problem. Scaling the momentum with a length 
scale £, we see that the instability signaled by the ap- 
pearance of a zero eigenvalue as t — > is accompanied by 
a diverging length scale 



£(* -► 0) oc (const, x t)- 1/4Vi 



(3.14a) 



Clearly, £ is the correlation length for the phase tran- 
sition that is the result of the instability of the normal 
metal phase. According to the usual definition of the 
correlation length exponent v, £ oc t~ v , one has 

v = 00 . (3.14b) 

We now turn to the order parameter susceptibility and 
the critical exponent 7. Since the diverging correlation 
length £ appears in the particle-particle spin-triplet chan- 
nel, we expect the order parameter to be one of the corre- 
sponding components of the Q matrix, e.g. \Q. The or- 
der parameter susceptibility will then be the correspond- 
ing two-poi nt cor relation function. In the symbolic nota- 
tion of Sec. [ill A , the susceptibility is given by 



X (VL n ) = (5 n |M-%) 



1,2,3,4 



<5l+2,n M 12 34 #3+4 in 



(3.15) 



Now con sider the eigenvalue equation for the matrix M, 
Eq. (3.2), and go into the basis of eigenfunctions /W to 
the eigenvalues X j. In serting two complete sets of eigen- 
functions in Eq. (|3.15]) we obtain 



(f^\8 n ) . (3.16) 



X (fi n )=^(5 n |/W) 



Let A be the smallest eigenvalue, and f\ the correspond- 
ing eigenfunction. Then the leading contribution to the 
diverging susceptibility as A — > is obtained by keep- 
ing only fx in Eq. (3.16). At zero external frequency we 
obtain 
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X (Q = 0) 



du> f\(uj) /A / dujfUuj) 



(3.17) 



The divergence of \ as A — > is determined by the ultra- 
violet behavior of the eigenfunction /. With Eq. ( 3.11 ) 
we obtain 



X(Q = 0) ocT 



(3.18a) 



This means that the critical exponent 7, defined by x cx 
i~ 7 , has its mean-field value, 



7 = 1 



(3.18b) 



Finally, we determine the critical wavenumber depen- 
dence of the order parameter susceptibility, i.e., the crit- 
ical exponent 77. Here the disorder part of the action 
make s a di fferenc e, and so we consider it explicitly. From 
Eqs. (3.4b) and (2.16a) we see that, with a = X> -1 , the 
eigenvalue A as a function of both t and the momentum 
p has the form A cx t 1 / 2 ^ + const, x p 2 . In a scaling 
sense,! we therefore have i 1 / 2 ^* 

2 



p 2 , or 



t 



In In 1 /p z 
lnl/p 2 



(3.19) 



plus terms that are less leading for p — > 0. Together with 
Eq. ( 3.18a ) this implies for the critical susceptibility 



x(P,t = 0) oc 



lnl/p 2 
In lnl/p 2 



(3.20a) 



According to the definition of the critical exponent 77, 
x(p,i = 0) cx p~^ 2 ~ v ^ this means 

T) = 2 , (3.20b) 

up to logarithmic corrections to power-law critical behav- 



IV. THE ORDERED PHASE 

A. The gap equation 

We now turn to t he or dered phase. Since the diverging 
susceptibility, Sec. Ill C| , is in the particle-particle spin 
triplet channel, we know that the order parameter will be 
a Q matrix in that sector. The symmetry group in that 
sector is U(l) x SU{2) (the usual gauge symmetry plus 
rotations in spin space), which is spontaneously broken 
to Z2 x U(l). We choose r = s = 1 as the direction in 
which ±hc symmetry is broken, and accordingly make an 
ansat£3 



r Q 12 (x) 



sp 



r-A- 12 ( x ) — <5ni,-ri 2 $rl #il A ni 



(4.1a) 
(4.1b) 



for the saddle-point values of the fields in the broken- 
symmetry phase. Using this ansatz in the saddle-point 
equations, Eqs. (2.9), yields the following equations for 
A n and Q n , 



1 



U3' 



+ ^-(A„) 2 



Here 



Kn 



V In 



n — m 



n + m 



and we have used the fact that A, 
function of the frequency. 



(4.2a) 



(4.2b) 



(4.2c) 



-A_„ is an odd 



It is clear from Eq. ( 2.7a ) (see the third term on the 
r.h.s of that equation) that A„ plays the role of the self 
energy in the particle-particle spin-triplet channel or of 
the gap function for the triplet superconductor. A crucial 
question now arises about the reality properties of A„. It 
follows from Eq. (4.2a) that A n is real (imaginary) if and 
only if Q n is real (imaginary). According to the formal 

( 2.6ej ), Q n should be 



2.6d) 



symmetry properties, Eqs 

real . However, as we have discussed at the end of Sec 
|lB|, this is misleading. In order for the field theory to be 
stable and correctly reproduce perturbation theory, the 
integration contour for the Q must be deformed, and for 
our saddle point to lie on the contour we must choose Q„, 
and hence A n , to be imaginary. Another requirement is 
of course that this choice minimizes the free energy. We 
will show in the next subsection that this is indeed-jthe 
case. We thus identify A„ = — iA„, with A„ real,ll3 as 
the gap function, and obtain the gap equation in the form 



V 



(4.3) 



Note that this gap equation has the same structure as 
the one in BCS theory (except for the frequency depen- 
dences of the kernel and the gap function), and that it 
isidentical to the gap equation that was derived in Ref. 
Ij.Ej It was shown in that reference that there are nonzero 
solutions for A„ for y > y c — 0.25, i.e. in the region in 
parameter space where we have found the normal metal 
to be unstable. We also note that the gap equation is the 
same irrespective of whether or not we explicitly keep the 
disorder term in the action. This is the triplet analog of 
Anderson's theorem. 

The gap equation, Eq. (4.3) was solved in Ref. 0, both 



analytically in certain limits and numerically for all fre- 
quencies, and there is no need to repeat this. Let us 
discuss, however, one important point that was not men- 
tioned in the earlier work. To this end, we first recall 
that by using the approximation , Eq . ( |3 .6| ) , for the ker- 
nel of the integral equation, Eq. (|4.3|), one can transform 
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FIG. 1. Low- frequency behavior of the gap function for 
y = 1.3 j/ c . See the text for an explanation of the numeri- 
cal procedure. 



the integral equation into an ODE, as we did for the 
eigenfunction in Sec. [II B. We briefly recall the most 
important features of the solution. On the imaginary 
axis, A„ = A(iujn) = S(L0 n ), and the limiting behavior of 
5(lo) for asymptotically small lu is 



5{u) — > 0) = — 2 yu> ln(w/w c ) 



(4.4a) 



with uj c a frequency scale that may depend on y, and for 
asymptotically large frequencies 

5{w — ► oo) = const, x cos (ji| 1/2 hxu + p\ , (4.4b) 

with a constant prefactor and a phase <f> that we could 
not determined analytically. Figure ^ shows the function 
5(u>) for small to, and Fig. g shows the oscillatory behav- 
ior at large a; on a logarithmic scale. For t his numerical 
solution at a fixed value of y, Eq. ( 4.4a ) was used to 
provide an initial condition for a very small but non-zero 
frequency, and then the ODE was iterated towards larger 
frequencies. For the initial condition, a value of lo c was 
chosen arbitrarily. It turns out that one finds a solution 
irrespective of the value of u> c chosen, and that all of these 
solutions are qualitatively the same. This behavior of 
the numerics is readily understood by inspecting the gap 
equation. It is easy to see that if 6(lo) is a solution, then 
so is 7(w) = b^ 1 6(bu>) with an arbitrary positive number 
b. This sc aling property holds for both the integral equa- 
tion, Eq. (4.3), and for the ODE that is derived from it 
by approximating the kernel. This ra ises t he question of 
what determines the scale u> c in Eq. (4.4a). In principle, 
the answer must be inherent in the gap equation itself, 
but we can obtain uj r from a much simpler argument. As 



we have seen in Sec. Ill , the transition is an ordinary con- 
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FIG. 2. Overall behavior of the gap function for y = 1.3 y c 
on a logarithmic frequency scale. See the text for an expla- 
nation of the numerical procedure. 



scaling will work. We further know that frequencies scale 
like wavenumbers squared, to ~ p 2 . From a scaling point 
of view, lo c is the critical frequency scale, and hence 



(4.5) 



where we have used Eq. (3.14a). This means that 
St^o (ui) —>■ for all u>, as one would expect, and the 
overall amplitude of the order parameter vanishes like 
u> c . For the critical exponent /3, defined by A cx t" , this 
implies 

P = 2v = oo . (4.6) 



B. The free energy 

We now check whether the saddle-point solution that 
we found to exist is stable, i.e. whether it has a free en- 
ergy that is lower than that of the normal metal state. 
In mean-field approximation, the free energy is given by 
the action, with the fields replaced by their mean-field 



values. By substituting Eqs. (4.1) in Eq. (2.7a) and ex- 
panding and resumming the Tr In term, we find for the 
free energy density in mean-field approximation 



2T y^^-t^Y. ln [i - W« + £)] . 



(4.7) 



where we have used Eqs. ( |4.2j ) to rewrite the term that 
was quadratic in Q n as a term that is bilinear in Q n 
and A n . We next use Eqs. (4.2) again to make this term 
quadratic in A n , 



tinuous phase transition with a diverging length scale, so 



7 



/ = 



T E^E ln [i-a^/k+^)] 



1 K t — ft 



(4.8a) 



where K 1 is the inverse of the operator K, Eq. (4.2c) 
in the sense that 



(4.8b) 



Now we use A„ = — iA n , and expand in powers of the 
order parameter. This yields the Landau expansion 



/ 



0(A 6 ) 



G„ m (k) = 6 (k) (r <8> s ) + S .F„(k) (n ® si) , 

(4.11a) 

where 

-(iw n + £k) 



^■„(k) = 



+ f, 2 + A 2 ' 

n ' ^>k ^ n 



-A„ 



(4.11b) 



ui+ai+Ai ■ 

Notice again the structural analogy to BCS theory. 
1. The anomalous susceptibility 



(4.11c) 



where 



t =^K~' 



Lx 1 \- 1 

T 6nm V^u J ? l +C ; 



and 



1 

2V 



E 



i 



> 



(4.9a) 



(4.9b) 



(4.9c) 



«„ is positive definite, and hence the theory is stable. 
Furthermore, t nm is negative for small K~ x or large K, 
i.e. for strong coupling. Hence the free energy in a state 
with a nonzero order parameter is indeed lower than in a 
state with zero order parameter, provided that the cou- 
pling is sufficiently strong, i.e. y > y c . All of these 
considerations are in exact analogy to the case of BCS 
theory. We note that if A n we re real, then the sign of 
the quadratic term in Eq. (4.9a) would be reversed, and 
consequently the free energy would favor order at small 
coupling, which is somewhat paradoxical. This, was actu- 
ally mentioned already in Berezinskiiis paper and later 
investigated in more detail by Heid.E3 



C. Linear response in the ordered phase 

We now calculate the linear response in the ordered 
phase, specifically the magnetic susceptibility and the 
electrical conductivity, in Gaussian approximation. This 
will establish that the ordered phase is really a supercon- 
ductor. 

To start with, we need the saddle-point Green function. 
It's inverse is (see Eq. (2.7a)) 



G- 1 = Gn 1 - iA 



(4.10) 



Using Eq. (4.1t), and inverting the 2x2 matrix, we ob- 
tain 



Before we calculate observable correlation functions, 
let us focus on the Goldstone modes. Although not di- 
rectly observable, they are of crucial impor tance for the 
former. As we already noted in Sec. IV A , the relevant 
symmetry group is U(l) x SU(2), which is spontaneously 
broken to the little group Z-i x U(l). The dimension of the 
relevant quotient space is dim (U(l) x SU(2)/U(l)) = 5, 
so there are five Goldstone modes. With our choice of the 
order parameter in the sector r = s = 1, they are repre- 
sented by the correlation functions 



^(;Q 1+nrl (k)jQ 2+Ilr2 (-k) 

1,2 



(4.12) 



with r — l;i — 1,2, or r = 2;i = 1,2,3. It will turn 
out that the r = 2, i = 1 Goldstone mode couples to 
the density response, so let us focus on the anomalous 
susceptibility 



Now consider the eigenvalue problem 

E ^12,34 /34 = A /12 . 
3,4 



(4.13) 



(4.14) 



This is the same eigenvalue problem as the one we con- 
sidered in Sec. HI in the symmetric phase. We therefore 
know that at zero momentum there is an eigenvalue Ao 
that is infinite everywhere i n the b roken symmetry phase. 
We further know (see Sec. Ill C| ) that the wavenumber 
dependence of A -1 quadratic. Perturbation theory can 
be used to convince oneself that the frequency depen- 
dence is also quadratic, although determining the pref- 
actors would be hard. Dropping all prefactors, we thus 
have the structure 



A( P ,n) 



i 



Q 2 



(4.15) 



everywhere in the ordered or broken symmetry phase. 
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2. The transverse current susceptibility 

Let us now calculate the transverse current suscepti- 
bility, which determines the magnetic susceptibility and 
hence the Meissner effect. By adding an appropriate 
source term to our action, we obtain for the transverse 
susceptibility in saddle-point approximation 

XT(k - 0) = -L J2 P 2 T£ H?„(p) 9n(p) 

p n 

+^ n (p)J-_ n (p)] , (4.16) 

in D = 2 di mensi ons. Notice that in a diagrammatic lan- 
guage, Eq. ( gig ) is a simple bubble. In the transverse 
channel, this is sufficient since all vertex corrections van- 
ish. In the language of the present field theory, saddle- 
point and Gaussian ap proxim ations coincide. Now we 
do the integrals in Eq. (4.16). A convenient way to do 



this is to put an upper cutoff on the integral over £ p , do 
that integral first, and let the cutoff go to infinity in the 
end. One can readily check that this procedure yields the 
correct answer in the BSC case (where convergence prob- 
lems can be avoided by doing the frequency integration 
first). We then findfe! 



XT(k-^O) 



duj 

» (ui 2 



3/2 



(4.17) 



The integral is positive definite, therefore xt < n/m and 
we have an ordinary Meissner effect. Notice that the re- 
ality properties of A and A are once again crucial for this 
conclusion: If A in the particle-particle spin-triplet chan- 
nel were real, then the correction to the free fermions 



result in Eq. (4.16) w ould be positive! We also note that 
5(lo) as given by Eq. (4.4a) renders the frequency integral 
in Eq. ( 4.16| ) infinite. This is an artifact of the theory that 
is due to a missing self-consistency. As was discussed in 
Ref. ||, taking into account the feedback of the nonzero or- 
der parameter on the electron-electron interaction would 
cure this problem. 



3. The longitudinal response 

The calculation of the longitudinal current, or density, 
response, and hence the conductivity, is much more in- 
volved. The reason is that in order to obtain a conserving 
approximation, one needs to solve the full Gaussian the- 
ory (in diagrammatic language, one needs to take the 
vertex corrections into account). This is well known in 
BCS theoryJiZI and the calculation for the present case 
proceeds analogously. As a check of our procedure, we 
have also performed the calculation for the BCS case and 
have reproduced the standard results. 



We start out by expanding the action, Eq. (2.7a) 
to Gaussian order about the saddle-point solution, Eqs 



(4.2). For the part that is quadratic in SA = A — A sp we 
find 



Tr (G sp SAG sp SA 



4 E ^EEE^w 



rs^l2,34 



(k)i(«5A) 34 (-k) 



(4.18a) 



Here 



«> 12i 3 4 (k) = 5 13 <5 24 ^(k)M r s % 

+5i 3 <5 2 ,-4^(k)M r s 1 ™ 01 



10 _io 



A-3Wi2(k)ACm( 
-*i 1 -4a 2 ,-3¥>8(k)M r 1 > 1 Mi 1 



(4.18b) 



with 4x4 matrices 



M w = 



, oi 



ZT 3 

-in 

-in o 
o in 

T 2 
-T2 



M 01 = 



M 11 = 



-in o 

o -in 

-ir 3 

-T 



in 



10 



T 2 

r 2 



(4.18c) 



and 



<P 00 V 
rnrn \ 



( k ) = -^E^(p)^™(p+ k ) . ( 4 - 18d ) 



and ip 01 , (p 10 , and ip 11 defined analogously with QQ in 
Eq. ( 4.18d| ) replaced by QJ-, TQ, and TT respectively. 

Similarly, the term that couples <5A and SQ can be 
written 



Tr (5A6Q) =4 J2 ^EEE^AMk) 



1,2,3,4 



k r,s i,j 



x %B (k)i(SQ) 3i (-k) , 



where 

ri^l2,3 4 ( k ) = S 13 S 24: S rs | 

The full Gaussian action is 
A G = Tr (g sp SAG sp SA) 4 Tr 



(sasq) 



(4.19a) 



(4.19b) 



(4.20) 



This Gaussian action is too complicated to be handled 
conveniently, and approximations are necessary. In or- 
der to make sensible approximations, we recall that our 
purpose is to calculate the density response in a conserv- 
ing approximation. In order to do so, it is crucial to 
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preserve the structure of the soft modes in the theory, 
while massive modes can be dealt with in very crude ap- 
proximations. It turns out that the soft mode structure 
is preserved by keeping only the imaginary parts of the 
fields in the r = 1,2,3 channels, and only the real parts 
in the r = channel, and integrating out the remaining 
parts in saddle-point approximation. This amounts to 
making the approximations 



;(<5A) 12 « \ 



12 



+ 
+ 



:{6A). 



1,-2 



ilm " 
vim 
ilm 
. Re , 



i(*A) 12 = l\ 



12 



(4.21a) 



and approximating 8Q by an analogously defined object 
q. A and q have the symmetry properties 



B12 



-— TtNf T(5i+2,3+4 ^12,34 



(4.22d) 



By adding an appropriate source term to the action 
it is easily checked that the expression for the density 
susceptibility \ in terms of the q is the same as in terms 
of the Q, viz. 

X (k,Q n ) = 16T^ ]T (°<Z 1+ „,i(k)°g 2+ri , 2 (-k; 

1,2 r=0,3 

(4.23a) 

Fro m th e expression of Q in terms of the fermion fields, 
Eq. (|2.2| ), it is easy to see that the contributions to Eq. 
( (4.23a|) from r = and r = 3 arc identical, except at zero 
external frequency. We can therefore write 



X{ 



(k, n n + 0) = 32T^ (§g 1+nil (k) ^ 2+n , 2 (-k) 



r' l l2 



-1,-2 ' 



(4.21b) 



(4.23b) 



+ ) r<?-l,-2 ' 



(4.21c) 



The net effect of these approximations is a symmetriza- 
tion of the action with respect to positive and negative 
frequency values. 

We now formally integrate out SA. This yields the 
action entirely in terms of 



and obtain the zero frequency susceptibility from Eq. 
( P~23b[) in the limit Q n -> Oil Combining Eqs. ( p~23b| ) 
and ( |4.22a| ) we have 



x(k,a i ^0) = 4T^™M i ;„ 1;2+ „ 2 (k) , (4.24a) 

1,2 

where 



M~ 



(4.24b) 



V 



p 1,2,3,4 r,s i.j 

\ij A^ 1 (r,\ _ ij U 
lrs^ 1 12,34VJ J i rs-^12,34 



iq u (-p) . (4.22a) 



Here l J s A 12 34 is the inverse of the matrix 



Instead of invertin g the m atrix M directly, it is conve- 
nient to rewrite Eq. (4.24b) as an integral equation, 



M _1 = A + ABAr 



(4.25) 



Ignoring the frequency indices for the time being, we need 
the matrix clement 



y a 

rs^ 1 12,34 



and 



+ 



+ , 
' + 



j 



y a 

rs^ 1 12,34 



y A 

rs-"--l, -2:3,4 



y 4 

rs^ 1 l,2;-3.-4 



y A 

rs^ L -l,-2;-3,-4 



(4.22b) 



33 



33^ + 32^-^23"^ 



(4.26) 



where we have used the structure of B, Eq. (4.22c) 
23M in turn obeys the integral equation 



23^ + 22^ B 23M 
10 a , 11 A a 10 



— 23^4 + 22^ {3 23A 

where (3 is the solution of 

(3 = B + B\\A(3 . 



(4.27a) 



(4.27b) 



( 3 is related to the anomalous susceptibility x^ a \ Eq. 
(PI, by 



US 



rs" 12,34 — u y 



with 



12.34 



(4.22c) 



(3 = HA \^B , 



(4.27c) 



with our approximation for that replaces Q in Eq. 
(EH by 9. 



10 



We are now in a position to calculate the density sus- 
ceptibil ity. D efining a vector \S n ) = #i,2-« we can write 
X, Eq. ( 4.24a| ), as a matrix element 



X (Kn n ^0) = 4T(S n \° 3 °M \k)\S n ) . 
Let us first consider the 'bubble contribution', 

xo(k,fin) = 4T(a n |§§ J 4(k)|<y n ) . 

Equations Q4.22b| ) and yield 

X0 (k,O„) = M% nA M - ^„,i( k )) 

' ,i(k) 



(4.28) 



(4.29a) 



VPi-. 



i(k))] 
(4.29b) 



Doing the integrals in the limit of small wavenumber and 
frequency yields 



Ji(k,n„) = jv F + o(k 2 ,n) 



(4.29c) 



Notice that the bubble contribution is not a conserving 
approximation for x, as the susceptibility does not vanish 
at zero wavenumber. As in the BCS case, particle number 
conservation or gauge invariance will be restored by the 
coupling of the density fluctuations to a soft or massless 
Goldstone mode, as we now proceed to show. 

We now consider the 'soft' part of x, i- e - the o n e tha t 
couples to the Goldstone mode. From Eqs. ( 4.28 ), ( 4.27 ) 
we find 

X S oft(k,0„) = (S n \ 3 lA(k)Bx^(k)BllA(k)\8 n ) 



5>„lMk)B|/) (f\x (a) (k)\f) 



x(f\Bl°A(k)\6 n ) , (4.30a) 



where we have inserted two complete sets {|/)} and 
{I/')} °f eigenfunctions of x ■ Since we are interested 
in the coupling to the soft mode, we need to keep only 
the eigenfunction / whose eigenvalu e A is infinite at zero 
momentum and frequency, see Eq. (4.15). We thus have 



Xs oft(k, Sl n ) = (6 n \° 3 lA(k)\g) A(k, Sl n ) {g\™A(k)\6 n ) , 

(4.30b) 

with |ff) = B\f) the transformed eigenfunction from 
Eq. (3.3). The structure of the eigenvalue we know, 
Eq. (4.15), and what remains to be done is to calcu- 
late the coupling integral. Using the symmetry property 
ffi2 = — ff2i of the eigenfunction that follows from Eq. 
(3.4a), we find 

(8 n \°lA(k)\g) = (g\l° 3 A(k)\5 n ) 

- [(fl+nA 9l+n,-l ~ Vi_„,i 9l-n-x) 
1 

+ (Pl+n,l 9l+n -1 ~ <Pi-n,l 9l-n,-l)] ■ 

(4.31a) 



h(k,Sl r 



Using the above symmetry property of g, and th e fact 
that the anomalous Green function J 7 , Eq. ( 1.11c ), is an 
antisymmetric function of the frequency, we see that I2 is 
also an odd function of its frequency argument. A Taylor 
expansion yields 



h(k,Sl r , 



-h(k, —Sl n ) oc fl r , 



(4.31b) 



We conclude that the density susceptibility at small fre- 
quencies and wavenumbers has the structure 



X (k,Sl n ) = N F 1- 



;Sll 



bk 2 



■ Sll 



(4.32a) 



where a, 6, c are constants that we have not determined 
explicitly. Since only the coupling to the Goldstone mode 
can restore the property x(k — > 0, Sl n ) = 0, we must have 
a = c. We have performed the analogous procedure for 
the BCS case, where a — c is easily confirmed explicitly. 
We conclude that the structure of the density suscepti- 
bility is the same as in the BCS case, viz. 



X(k,0 r , 



Nv 



{vk) 2 



sil + W 2 



(4.32b) 



with v the velocity of the sound-like mode that is 
the analog of the, Anderson-Bogolubov mode in BCS 
superconductors.cij Here we have not determined v ex- 
plicitly. 

The electrical conductivity a is determined by x y i a 



a(Sl n ) = e 2 lim — £ x(k, Sl n ) 
k-»0 k z 



(4.33) 



Combining Eqs. Q4.32b| ) and ( |4.33| ), we see that a(Sl n ) oc 
1/Sl n , and the real part of the conductivity as a function 
of real frequencies has a delta-function contribution 



e 2 N F v 2 
Re a (Sl) = —S(Sl) 



(4.34) 



The ordered phase is thus a true superconductor. If we 
neglect the explicit dis order term in the action, then the 
delta function in Eq. ( 4.34 ) will exhaust the f-sum rule, 
which in turn determines the velocity v. However, with 
the disorder term taken into account, there is a continu- 
ous spectrum in addition to the delta function, and the 
determination of v would require an explicit calculation 
of the prefactors in Eq. (4.32a). 



V. DISCUSSION AND CONCLUSION 

A. Summary of results 

Let us summarize. We have considered in detail the 
consequences of the attractive interaction in the particle- 
particle spin-triplet channel that has been predicted to 
exist in 2-D disordered interacting electron systemsJj We 
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have found that this interaction leads to an instability of 
any normal metal phase in the system. The instability is 
characterized by a diverging length scale £, and has all 
the characteristics of a continuous quantum phase tran- 
sition. The critical exponents v 1 7, and 77 are 



7 = 1 



V = 2 



(5.1a) 



where the infinite value of v indicates an essential singu- 
larity in the dependence of £ on the distance from criti- 
cality t, 



£cx (ct)- 1 / 4 ^ 



(5.1b) 



with c a constant, and the value of r\ is to be understood 
as a complicated logarithmic dependence of the critical 
order p arameter susceptibility on the wavenumber, Eq. 
(|3.20a|) . 

This instability identifies the order parameter of the 
phase transition, viz. an anomalous density in the 
particle-particle spin-triplet channel. Approaching the 
transition from the ordered phase, the order parameter 
critical exponent is 



(3 = 2v 



(5.2a) 



in the sense that the amplitude of the frequency depen- 
dent order parameter vanishes like 



a oc r 



(5.2b) 



The phase transition thus has all the characteristics of 
an ordinary quantum critical point. The ordered phase 
has all the characteristics of a superconductor, including 
a Meissner effect and an infinite static electrical conduc- 
tivity. While these conclusions are general, the detailed 
critical behavior, and in particular the values of the expo- 
nents, depe nd on our particular choice of the interaction 
kernel, Eq. (|2~8^ ). 



B. Discussion 

As we have mentioned before, the reality properties 
of the gap function are crucial for the stability of our 
superconducting state. If the expectation value of \Q 
were real, then the sign of the (A m ) 2 term in Eq. ( |4.3| ) 
would be negative. Irrespective of whether or not this 
modified gap equation has physical solutions, this would 
lead to a free energy that is minimized by a non-zero 
gap function at,weak coupling. This was noticed already 
by Berezinskii,L!l who proposed that odd-gap supercon- 
ductivity or superfluidity would appear in a temperature 
window. Later, a_more detailed discussion of the free 
energy was givenjl3 and it was shown that the different 
sign in the gap equation_also leads to the Meissner kernel 
having the wrong sign.Ej This led to n the abandonment 
of Berezinskii-type odd-gap proposalsB in the context of 
high-T c superconductivity. Instead, these authors then 



considered more complicated scenarios of composite or- 
der parameters that do not suffer from these problems. 

In the present theory, these problems are solved by re- 
alizing that there is a choice with respect to the reality 
properties of the gap function we call A„, and that the 
physical choice is to make it imaginary. Let us recall 
the observations that led to this realization. First of all, 
with the conventional choice of a real A n , the field theory 
that results from expanding about a Fermi liquid sad- 
dle point is unstable in the particle-particle spin-triplet 
channel, which forces a deformation of the integration 
contour. The validity of this procedure was then checked 
by comparing the perturbation theory that results from 
the field theory with ordinary perturbation theory within 
the framework of a second quantized Hamiltonian. The 
comparison can be made in the very simple case of non- 
interacting electrons, where there is no question as to 
what the correct answer is. The perturbative result is 
also in quantitative agreement with experiment. Next, 
in order for the saddle point that describes the supercon- 
ducting phase to lie on the deformed contour, the saddle- 
point fields, which correspond to anomalous expectation 
values in a fermionic formulation, must be chosen to be 
imaginary. This choice turns out to minimize the saddle- 
point free energy, while the conventional choice would 
maximize it. The conclusion is that our interpretation of 
(\Q n -i - n ) as imaginary is correct, and there is no ambi- 
guity or freedom left. We note that the particle-particle 
spin-triplet channel is the only one where the determina- 
tion of the reality properties of the expectation values is 
not straightforward.Eil In all other channels, the symme- 
try properties of the fluctuating fields, Eqs. (2J3), agree 
with the convergence requirements of the field theory, 
and no deformation of the contour is necessary. We do 
not know why the particle-particle spin-triplet channel is 
special in this respect. 

One might object tha t the symmetry properties of 
the Q matrix, Eqs. d2.6| ), just encode the properties of 
the underlying fermionic fields, and must therefore also 
hold after taking expectation values, which would make 
(iQ n -i _„) reai - This objection is fallacious, however, 
due to the following reasons. The symmetry properties 
hold for the fluctuating fields, which are just dummy in- 
tegration variables and have no direct physical meaning. 
Deforming the integration contour to apply the method 
of steepest descent is therefore permissible as long as one 
does not cross any singularities. This we have checked be 
comparing with fermionic perturbation theory. The ex- 
pectation value (\Q n _ 1 ), on the other hand, does have 
a physical meaning, but it does not follow from Eqs. ( |2.6| ) 
that it must be real. In the absence of spontaneous sym- 



metry breaking, of course (\Q r _ 



0. To obtain 



a non-zero expectation value, one must consider a small 
field that is conjugate to \Q n _ 1 . If this field were 

static, then hermiticity would indeed require (\Q n _ 1 __) 
to be real. However, in the particle-particle spin-triplet 
channel the conjugate field cannot be static. The real- 



12 



ity properties of the anomalous expectation value then 
depend on the properties of the unphysical, time depen- 
dent, symmetry breaking field, which are not given a pri- 
ori. They could be determined from the requirement that 
(\Q n -i - n ) be imaginary, but we have not done that. We 
also note that our sad dle po int violates only the hermitic- 
ity r equir ement, Eq. (2.6e), but not the Pauli principle, 
Eq. (2.6d), in agreement with the above discussion. 

One might wonder why the static conductivity is af- 
fected by the existence of a gap function, given that the 
latter is zero at zero frequency. The answer lies in the fact 
that the infinite conductivity, like the Meissner effect, is 
a direct consequence of the broken symmetry and the re- 
sulting Goldstone modeB3 Indeed, in BCS theory with its 
constant gap function, the prefactor of the delta-function 
in the conductivity is independent of the gap. More 
specifically, given the existence of a sound-like Goldstone 
mode that couples to the density fluctuations, particl e 
number conservation determines the form, Eq. (4.32b), 
of the density susceptibility, which in turn guarantees an 
infinite static conductivity. On an even more technical 
level, we o bserve that while the Goldstone susceptibility, 
Eq. (4.12), indeed vanishes at zero external frequency, 
the actual coupling to t he de nsity fluctuations is more 
complicated. From Eq. ( 4.31a ) we see that the coupling 
i s giv en by the more general susceptibility Xi2 34j Eq. 
( 4.13 ), multiplied by the cigenfunction g, which is a gen- 
eralized gap function. The net result, as far as frequency 
dependences is concerned, is the same as for conventional 
superconductors. 

We have neglected the explicit disorder term in the ac- 
tion, and have taken the quenched disorder into account 
only implicitly, by using the particle-particle spin-triplet 
interaction from Ref. f|. At the level of the Gaussian 
theory, keeping the disorder explicitly just complicates 
the calculations without leading to any essential changes, 
as we have pointed out throughout the paper. For in- 
stance, the gap equation is independent of the disorder 
(this is the spin-triplet analog of Anderson's theorem), 
and the critical behavior is qualitatively unchanged. An- 
other (minor) modification is that the (5-function contri- 
bution to the conductivity no longer exhausts the f-sum 
rule, but rather sits on top of a continuous spectrum. 
This conclusion from the Gaussian theory would indeed 
be correct if w e could work in dimensions D > 2. How- 
ever, Eqs. ( 2.8 ) are valid only in D = 2, and in D > 2 the 



analogous interaction does not lead to a superconducting 
ground state.Q In D — 2, on the other hand, the explicit 
disorder term in the action has more subtle, qualitative 
consequences, which we discuss next. 



C. Outlook 

Throughout this paper, we have ignored the well- 
known fact that, in a 2-D disordered interacting elec- 
tron system at zero temperature, the paramagnetic nor- 



mal metal phase is also unstable against the formation 
of an insulator, and a ferromagnetic phase. As a result, 
as we have stressed elsewhere, E2l the nature of the ground 
state in the absence of spin-flip scattering processes is not 
known, even in the absence of the attractive spin-triplet 
interaction discussed above. 

We conclude that in a 2-D disordered electron fluid, 
there are at least three competing instabilities of the 
normal metal ground state. Which one of these wins 
is a priori unclear, and answering this question would re- 
quire a detailed and consistent study of the free energies 
for the various phases, which has not been attempted so 
far. It is certainly conceivable that different instabilities 
dominate in different parameter regimes, leading to very 
different ground states. It is also possible that the ten- 
dencies towards forming an insulator and a spin-triplet 
even-parity superconductor, respectively, effectively can- 
cel one another, leading to a normal metal-like phase at 
least in a certain temperature window. This last possi- 
bility is the basic idea behind recent proposals to explain 
the observed metal-insulator transition in Si MOSFETs 
and other 2-D electron systems in terms of our exotic 
superconductivity.B However, a better understanding of 
these issues will require much more work. 
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APPENDIX A: WEAK-LOCALIZATION 
CORRECTIONS, AND THE STABILITY OF THE 
FIELD THEORY 

The purpose of this appendix is to check that our field 
theory, defined by th e act ion and the integration con- 
tour specified in Sec. II B, correctly reproduces the re- 
sults of standard perturbation theory. For this purpose, 
we choose to calculate the well-known weak-localization 
corrections to the conductivity within our formalism. 

To this end, we restore the disorder part of the action, 
^4dis, although we will not need it explicitly. It is well 
known that non-interacting disordered electrons in D = 2 
dimensions in the absence of spin-orbit scattering do not 
have a metallic phase.C3 In perturbation theory, start- 
ing from the Boltzmann value for the conductivity a and 
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performing a disorder expansion, this phenomenon man- 
ifests itself as corrections at first order that diverge log- 
arithmically as the frequency fl approaches zero. These 
corrections are of the form 



5a = (a s + at) In Qt 

e F r 



(Al) 



Here ep is the Fermi energy, r is the elastic scattering 
mean-free time, and the coefficients a s and Ot are due to 
contributions from the particle-particle spin-singlet and 
spin-triplet channels, respectively. We need only the ra- 
tio of these two coefficients. In the absence of any spin- 
flip scattering processes as well aSr-ef magnetic fields, 
fermionic perturbation theory yieldsc3 



a t = — 3a s 



(A2) 



This result has,been confirmed by experiments on thin 
metallic films.Ej By deliberately contaminating the sam- 
ple with a strong spin-orbit scatterer, the triplet contri- 
bution to the weak-localization effect can be eliminated, 
which leads to a sign change of the effect as well as a 
reduction of the amplitude by a factor of 1/2. Eq. ( [A2] ) 
is thus very well confirmed. _ 

Perturbation theory within the Q-field theoryQ yields 
structurally 

Saoc [(-°W)(P)?W)(-P)> trsosotrrm 



-3(\(6Q)(p)\(6Q)(-p)) trs^tmn] 



^E[(iW)(p)?w)(-p)) 

-3(l(5Q)(p)l(SQ)(-p))] 



(A3) 



For notational simplicity we have suppressed the fre- 
quency labels on the Q matrices, which are related to 
the external frequency f2. In the presence _pf disorder, 
the (Q Q) correlation functions are diffusive,0 and there- 
fore the moment um integr al y ields a In fir in D = 2. 
Comparing Eqs. (Al) and (A3) we see that the field the- 
ory yields the same result as fermionic perturbation the- 
ory, provided that the two correlation functions have the 
same sign and magnitude. This is in deed the case with 
the procedure explained in Sec. II B, or, equivalently, if 
one ignores any convergence questions and formally per- 
forms the Gaussian integrals over Q, as was done in Ref. 
^. We conclude that the field theory, with the particle- 
particle spin-triplet matrix elements integrated over the 
contour chosen in Sec. II B, correctly reproduces pertur- 



bation theory, which in turn is in quantitative agreement 
with experiment. 
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